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Abstract 

We consider the expulsion of the magnetic field from the super-conducting core 
of a neutron star and its subsequent decay in the crust. Particular attention is 
paid to a strong feedback of the distortion of magnetic field lines in the crust on 
the expulsion of the flux from the core. This causes a considerable delay of the 
core fiux expulsion if the initial field strength is larger than 10^^ G. It is shown 
that the hypothesis on the magnetic field expulsion induced by the neutron star 
spin-down is adequate only for a relatively weak initial magnetic field B ~ 10^^ G. 
The expulsion time-scale depends not only on the conductivity of the crust, but 
also on the initial magnetic field strength itself. Our model of the field evolution 
naturally explains the existence of the residual magnetic field of neutron stars. Its 
strength is correlated with the impurity concentration in neutron star crusts and 
anti-correlated with the initial field strengths. 
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1 Introduction 

Where the neutron star magnetic field (MF) is located, by which processes it is generated 
and what determines its strength, structure and evolution, is the subject of scientific 
debates since neutron stars are conceivable for human mind. Regarding the field location 
there exist basically two qualitatively different ideas: the field is either present in the 
entire star, that is, in the core and the crust, or is located mainly in the crust. Whether 
the entire neutron star or only the crust is penetrated by the field depends strongly on 
the mechanism responsible for its generation. 

Both the simple idea of fiux conservation during the gravitational collapse and the 
action of a dynamo in the convective proto-neutron star (Thompson & Duncan 1993) 
result in a field structure which penetrates the entire star. On the other hand it is 
possible that fall-back accretion after the supernova submerges any initially existent field 
(see Geppert, Page & Zannias 1999). Since the re-diffusion process may last a long time, 
the existence of young pulsars (PSRs) could be explained - at least to some extent - 
by the action of a thermoelectric instability (Blandford, Applegate & Hernquist 1983, 
Urpin, Levshakov & Yakovlev 1986, Wiebicke & Geppert 1996). A strong MF can be 
produced rather rapidly in the surface layer of the neutron star by transforming heat flux 
into magnetic flux. If so, the surface MF which governs e.g. the spin-down of the PSR 
is maintained by currents conflned in the crust alone. However, the currently existing 
models describe only the ampliflcation of large scale toroidal fields, the production of the 
observed dipolar MF is still an open problem. 

Although the assumption of a purely crustal MF explains well the observed long- 
term evolution of isolated neutron stars (Urpin & Konenkov 1997) and of neutron stars 
accreting in binary systems (Konar & Bhattacharya 1997, Urpin, Geppert & Konenkov 
1998, Urpin, Konenkov & Geppert 1998) the evolution and effect of a core MF have to be 
investigated too. The existence of a core MF may explain transient features like the post- 
glitch behaviour (Alpar, Langer & Sauls 1984, Chau, Cheng & Ding, 1992) as well as the 
presence of residual surface flelds in old neutron stars (e.g. Jahan Miri & Bhattacharya 
1994). 

It is generally accepted that a phase transition into the super-fluid/super-conductive 
state takes place in the core rather early in a neutron star's life if the temperature drops 
below 10^° K (see Alpar 1991, Page 1998). Given the typical neutron star MF strength, 
the core behaves as a type II superconductor (Baym, Pethick & Pines, 1969), i.e. the MF 
penetrates the core in quantized flux tubes (fluxoids). The evolution of the core fleld is 
connected with movements of the fluxoids caused by forces acting upon them. Under the 
action of these forces the fluxoids move towards the core-crust boundary, the magnetic 
flux is expelled from the core into the crust, where the magnetic field suffers Ohmic decay. 

Ding, Cheng & Chau 1993 (hereafter DCC) described these forces and developed a 
method to calculate the radial velocity of the fluxoids, which in turn determines the rate 
of flux expulsion from the core. There are the buoyancy force (Muslimov & Tsygan 1985), 
the drag force (Harvey, Ruderman & Shaham, 1986) and the force exerted by the neutron 



vortices which act upon the fluxoids. The interaction between vortices and fluxoids is 
determined by their pinning energy; in the course of neutron star spin-down the vortices 
move outward and fluxoids can move faster, comove or move slower than vortices. 

However, the motion of the proton flux tubes in the core leads to a distortion of the field 
structure in the crust near the crust-core boundary. This changes the magnetic energy 
in the crust which, in turn, is the source of a force that influences the flux expulsion. 
The aim of this paper is to consider that "back reaction" of the crust onto the fluxoid 
movement. 

The paper is organized as follows. In section 2 we describe the forces acting on the 
fluxoids in the core, the procedure of self-consistent determination of the inner boundary 
condition for the crustal field decay, and how to find the surface and core MF as func- 
tions of time. Section 3 represents the numerical results and in section 4 we discuss the 
consequences of crustal effects for the long term evolution of isolated neutron stars and 
compare our results with those of other authors. 

2 Description of the model 

As shown by Alpar, Langer & Sauls (1984) (see also DCC) both the proton fluxoids and 
the neutron vortices are associated with large magnetic fields of comparable strengths. 
Interaction of both structures involves a pinning energy of Ep ~ 10 MeV per intersection, 
and the strongest pinning occurs in the core regions where the flow of vortices is nearly 
perpendicular to the fluxoids. 

The flow velocity of the vortices, Vn, is governed by the spin-down torque, which in a 
steady state is assumed to be determined by the magnetic dipole braking of the neutron 
star (see DCC): 

Mt) - «« , (1) 

where 

m = 5^M« ^ (2) 

Here, Qs is the angular velocity of the core super-fluid. Its difference to the observed 
rotation rate of the crust, Qc, causes the driving Magnus force; R is the radius of the 
neutron star, x is the angle between the rotational and the magnetic axis of the neutron 
star, / its moment of inertia and B^ is the surface magnetic field in the equatorial plane. 
We will assume that equation (1) holds for all distances from the rotation axis r < R^ 
where Re is the radius of the super-conducting core. 

With respect to the relation between the radial velocity of the proton fluxoids, Vp, 
(for justifications to neglect their azimuthal motion see discussion in DCC) and f„ one 
can distinguish three different regimes. Either the vortices cross the fluxoids {forward 



creeping), or both types of tubes move with the same velocity {comoving) , or the fluxoids 
move even faster than the vortices {reverse creeping). Thus, in case of forward creeping 
the vortices exert a force onto the fluxoids which drives them outward, while in case 
of reverse creeping a resistive force acts upon the fluxoids which counteracts the flux 
expulsion. It is clear that this vortex acting force /„ depends on the angular velocity lag 
uj = Qs — ^c and changes its direction, when the Magnus force changes its sign because 
the core super-fluid rotates slower than the crust and the charged components of the core 
coupled with it. The maximum lag that can be sustained by the pinning force, Ucr oc Ep, 
defines also the maximum force a vortex can exert onto a fluxoid. This force per unit 
length of the fluxoid is given by (DCC, Chau, Cheng & Ding, 1993) 

_ 2%prQs{t)uj{t) 

where $o denotes the quantized flux per fluxoid, $o = hc/2e = 2 ■ 10^^ G-cm^, p is the 
density of the matter and B^ is the mean strength of the core MF, proportional to the 
number of fluxoids in the core Np and given by Be = ^oNp/nRl. 

Independently of the rotational evolution of the neutron star the fluxoids are affected 
by a buoyancy force ft,, which is caused by the magnetic stress at the surface of the fluxoid. 
This force per unit length of the fluxoids is given by (Muslimov & Tsygan, 1985): 



The relation between the London penetration length A and the coherence length C, deflnes 
the condition for the formation of a type II superconductor, it requires that ^ < \/\/2. 

The drag force is caused by the scattering of the degenerate ultra-relativistic electrons 
at the fluxoid MF. In an approximation which neglects collective effects of the fluxoids, 
the drag force per unit length is proportional to Vp and is given by (Harvey, Ruderman & 
Shaham, 1986): 

_ Sn uee'^^l vp , , 

64 EpX c 
where Ue is the number density of the electrons in the core, about 5% of the neutron 
number density, and Ep is the Fermi energy of the electrons. For the Fermi energy as well 
as for A and ^ we take the values determined by the density at the crust-core interface, 
Pc- Collective effects which modify the drag force become important only for Be > 10^^ 
G, when the mean fleld becomes comparable to the fleld inside the fluxoids. 
Additionally, DCC take into account tension forces which result in a coefficient of fb of 
the order of unity in almost the entire star. For simplification we will assume that the 
fluxoids are not bended but parallel to each other while moving together with their roots. 
According to DCC, the velocity of the fluxoids is derived from fn + fb + fv{vp) = 0. 
However, the MF which is concentrated in flux tubes in the core, penetrates the crust 
too. The movement of the roots of the fluxoids leads to a bending of the magnetic fleld 



lines in the crust (see figure 1). The power done by the core forces is equal to the net 
Poynting flux through the core surface, which, in turn, is equal to the rate of change of 
the sum of mechanical and field energy in the volume outside the core. We will assume 
the crust to be crystallized so that the change of mechanical energy corresponds only to 
the Joule heating of the crust. The contribution of the electric field to the field energy 
can certainly be neglected. Thus, we find for the balance of powers: 

E [ifn + fb + f.)v,dl= nl<W + ^ l^dV . (6) 

nuxoids T/ T/ 

Vcrust y 

Note, that f^ depends on Vp as well as the r.h.s. of equation (6). In the l.h.s. the 
integration is performed over a fluxoid's length, the summation runs over the number 
of fluxoids. For sake of simplicity we consider all terms of the l.h.s. of equation (6) to 
be position independent, take them at r = i?c and approximate that side by (/„ + /t + 



fv)Np (1) Vp. Since the length of a fluxoid is / = 2jRl — r^, its mean value is (/) = ARc/3. 
The number of fluxoids decreases with core field decay. The multiplication with the fluxoid 
velocity Vp yields the core forces power which has to be balanced by Joule heating and 
the change of the magnetic energy. While the Joule heating is restricted to the volume of 
the crust, the change of the magnetic energy has to be calculated both for the crust and 
the vacuum environment of the neutron star. We introduce the core forces 



n,b,v 



U,b,. ■ ^Rc/S ■ Np (7) 



and the crustal force 



crust 



^ldV + ^ f — dV\. (8) 

v„\ J a dt i Stt ' ^ ^ 



Thus, equation (6) can be rewritten in the form F„ + F^ + F„ + Fcmst = 0. 
In this paper we consider the evolution of a purely poloidal dipolar MF. Hence, the 
appropriate representation for the MF is that of the Stokes stream function S{r, t) which 
is related to the vector potential A = (0, 0, A^) with A^ = S{r,t) sm9/r, where r and 9 
are the spherical radius and polar angle, respectively. Then, we can express the spherical 
field components, Br and Bg, in terms of S{r,t), 

„ 2,5 ^ „ sin^ OS 

B, = —cos9, Be = — . 9 

We normalize S{r, t) by its initial surface value at the equator, S{R^ 0) = BeoR"^, s{r, t) = 
S{r, t)/S{R, 0), so that Be{t) = B^q ■ s{R, t). 

While in the core there is no ohmic decay but motion of the fluxoids, ohmic diffusion 
determines the field evolution in the solid crust. Thus, the field evolution in the core 
(r < Re) is governed by 



but in the crust {Re < r < R) hj 

ds _ c2 /dh _ 2s\ 
dt 47ra ydr"^ r'^ J 

Therefore, we find for the r.h.s. of equation (6) the following expressions: 
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dt J 8tt ^° dt . ^ . , 

Vcrust \ -Re R, 
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dt J 8n 3 ^ ' dt ^ ' 

V vacuum 

The assumption of a homogeneous mean MF in the core leads to the following expression 
for the stream function there 

s{r,t) = ^r^ (15) 

The electric conductivity in the solid crust consists mainly of contributions from electron- 
phonon and electron-impurity scattering. Electron-phonon interactions dominate the 
transport at high temperatures and relatively low densities, whereas the impurity con- 
centration determines the conductivity at lower temperatures and larger densities. We 
use the numerical data for the phonon conductivity obtained by Itoh, Hayashi & Koyama 
(1993) and an analytical expression for the impurity conductivity derived by Yakovlev & 
Urpin (1980). The crustal temperature which infiuences the electron-phonon conductiv- 
ity is taken from cooling curves calculated for different neutron star models by Van Riper 
(1991). For the chemical composition we adopt that of cold catalyzed matter. Note that 
the longterm crustal field decay in isolated neutron stars is determined by the impurity 
concentration mainly which is characterized by the impurity parameter Q because after 
~ 10® years the neutron star cools down completely. Moreover, the currents generated 
by fiux expulsion from the core in the crust are located in its deep high density regions. 
Thus, this field evolution is almost insensitive to the cooling history of the neutron star. 
The assumed homogeneity of the core MF leads to following ansatz for the fiuxoid velocity: 

Vp = a(t)r. (16) 

Thus, the solution of equation (10) is 



s{r, t) = s{r, 0) exp (-2a(t)t), (17) 

where a{t) has to be determined by the self-consistent solution of equation (6). In that 
way, Vp determines the inner boundary condition for equation (11). 
Since the MF in the center of the neutron star has to be regular we demand as the inner 
boundary condition that s{r,t)/r'^ remains finite with r ^ 0, which is fulfilled by our 
choice of s{r,t) in the core. The boundary condition at the neutron star surface has to 
ensure a matching of the interior field with the external dipolar field. That is expressed 
by the requirement i?|^ = —s at r = R. 

When solving equation (11) the boundary condition at the crust-core boundary is given 
by equation (17) with s{Rc,t), which decreases generally non-exponentially because a 
depends on time. 

In order to solve equation (6) we follow the ideas of DCC, however, we have to consider 
additionally the crustal effects given by the r.h.s. of the equation (6). By the choice of 
the strength of the core MF, a profile of s{r, 0) in the crust, an initial rotational period 
and a certain EOS which determines mass, radius and moment of inertia of the neutron 
star, we define the values of the parameters at t = 0. 

Then we start with Vp = f„ or a{t) = k{t)Vt^/2 where k{t) is given by equation (2). 
Adopting that equality we calculate s{Rc,t) from equation (17) which represents the 
inner boundary condition for equation (11). By use of equations (12), (13) and (14), 
we calculate the r.h.s. of equation (6), or Fcmst- After calculating Ff, and F^,, finally, 
we obtain F„ and thus uj. The core MF defines also the maximum lag ujcr that can be 
sustained by the pinning force (see equation (11) in DCC), which gives an upper limit 
of Fn- If —i^cr < UJ < uJcr then the fiuxoids and vortices comove and Vp, found by the 
procedure described above, is the real velocity of the fiuxoids. 

li u > UJcr then the fiuxoids are in the forward creeping regime. Since Ucr defines the 
maximum force which can be exerted onto the fluxoid we set uj = uJcr and calculate 
FniuJcr)- The solution of equation (6) with FniuJcr) provides Vp. At the forward creeping 
stage Vp < Vn, thus Vp will be found as a root of equation (6) in the interval [0, f „] using the 
bisectional method, li uj < —uJcr then the reverse creeping mode is reached. In this case 
we set u = —UJcr, calculate F„(— cUcr) and solve again equation (6). This yields Vp > f„. 
After having determined the velocity of the fiuxoids in that way we use equations (11) 
and (17) to calculate the evolution of s(r, t) and to obtain the temporal behaviour of the 
MF at the crust-core boundary Be and at the neutron star surface B^. From the latter we 
deduce the evolution of the spin period P{t) = 27r/r2. The described set of computations 
is performed at each time step. 

3 Numerical results 

Although the induction equation is linear, the problem is a nonlinear one in terms of the 
magnetic field. Additionally, equation (6) which determines the flux expulsion velocity 
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is nonlinear in Vp too. Thus, we have to consider the time evolution for different sets of 
initial values of B^ and B^ and to investigate their temporal evolution. 

For the calculations below we use the standard cooling scenario considered by Van 
Riper 1991 for the neutron star model based on the Friedman-Pandharipande (FP) equa- 
tion of state (Friedman & Pandharipande 1981) with the total mass M = IAMq, radius 
of the star R = 10.61 km and radius of the core Re = 9.67 km. 

In the present paper we consider the case that the homogeneous core MF has initially 
the same strength as the surface MF (see curve 1 in figures 1, 2). At first we assume 
for both Bf. and B^ the standard values 10^^ G and for the initial spin period 10 ms. 
As described above, the expulsion of the core MF causes distortions of the crustal field 
in the vicinity of the core-crust boundary (see figures 1, 2), generating currents just in 
that region. Hence, the decay of the crustal field, characterized by the decay time of 
the surface MF r^, is almost completely determined by the impurity parameter Q. The 
larger Q the more rapidly the crustal field may decay which, in turn, makes the braking of 
neutron star rotation by magneto-dipole radiation less efficient. If the MF of the neutron 
star is confined in the crust only Tg ~ 10^ Q~^ years is the appropriate estimate for the 
chosen neutron star model (Urpin & Konenkov 1997). However, the existence of a core 
MF and its expulsion into the crust causes a nonlinear dependence of the inner boundary 
condition for equation (11) on the MF evolution and deviations from the above estimate 
for Tg will occur. 

In figure 3 we show from top to bottom the evolution of the MF, of the velocities of both 
types of flux tubes, of the forces acting upon the fluxoids, and of the rotational period of 
the neutron star, calculated by the procedure described in the section above. Calculations 
were performed for three values of the impurity parameter Q = 1 (left column), Q = 0.1 
(middle column) and Q = 0.01 (right column). 

Almost independently on Q, the velocity of the fluxoids at times t < lO'^'^ is determined 
by the balance of the vortex acting force -F„ and the drag force F^. Since the spin-down of 
the neutron star is very effective, the vortices move outward fastly thereby cutting through 
the fluxoids: they are in the forward creeping regime. At this stage Vp falls down from 
10~^ to about 10~* cm/s. However, the characteristic time of expulsion Tg = Rc/vp ~ 10® 
years, and the core MF remains almost constant. Nevertheless, even during that relatively 
short period a strong meridional fleld component Bg is generated in the crust close to the 
crust-core interface. This gives raise of the crustal force Fcmst, and later Vp is determined 
by the balance of crustal and buoyancy forces. The power-like decay of -F„ simply reflects 
the increase of the spin period of the neutron star P oc \/i with an almost constant 
surface MF and the subsequent decrease of the vortex number density in the core. At 
this stage, Vp is dependent on Q: the lower Q the lower Vp, an effect of the crustal force, 
counteracting the expulsion. As long as the fluxoids are in the forward creeping regime. 
Be remains nearly constant. The forward creeping regime is followed by the comoving 
regime during which F„ changes its sign, and vortices and fluxoids move with the same 
velocity. Until that moment, the flux expulsion rate is governed by the balance of crustal 



and buoyant forces, however, the crustal force becomes ineffective on larger time-scales, 
t > Tg. This is the reason for the increase of Vp and, hence, the decrease of B^. Note, 
that the evolution of Be follows closely that of Be- The expulsion lasts till another force 
can prevent it. It appears to be a common situation that at late evolutionary stages 
the velocity of the fluxoids is determined by the balance of buoyancy and vortex acting 
forces. The amount of flux which has to be expelled from the core in order to reach such 
a balance is easy to predict by comparison of F„ and F), just before the beginning of the 
expulsion. Since Fb oc Be but Fn oc y/B^, the residual MF strength can be estimated by 




BeO \Fb{Te 

The longer the time-scale of expulsion (due to the resistive effect of the crustal force), 
the lower is the residual field strength. This is because the number of vortices and Fn 
are lower if the neutron star spins down to a greater rotational period during a longer 
time. Thus, -F„ will be able to balance Fi, on a lower MF level. In turn, the time-scale 
of expulsion depends on the impurity concentration. Namely, for Q = 1, 0.1, 0.01 the 
residual field strength is about 10^°, 10^, 10^ G, while Tg is about 10^, 10^, 10^ years, 
respectively. The velocity of the fluxoids increases just after the beginning of the reverse 
creeping stage, and decreases sharply when Ff, balances F„. Concluding, we want to point 
out that the main force which expels the fluxoids from the core is the buoyancy force. 
The vortex acting force influences the evolution of the MF of very young and very old 
neutron stars, while the effect of the drag force is restricted to the very early stage only. 
In figure 4 we show the same quantities as in figure 3 but the initial strength of the core 
and surface MF is assumed to be only 10^^ G. This leads to a drastically less efficient 
spin-down of the neutron star. Because the magnetic energy stored in the crust and in 
the surrounding vacuum is oc B'^, the influence of the crustal forces on the expulsion is 
much weaker than in the case of -Beo = 10^^ G. For a comparatively low fleld strength 
the situation is similar to that considered by DCC because practically only Fb, Fn and 
-Fj, determine the fleld evolution. Due to the much slower spin-down, from the beginning 
both vortices and fluxoids are moved with the same velocity. Again, in the young neutron 
star (t < 10^ years), Fn and F^ are much greater than Fcmst and Fb, while in the old 
neutron star Fn and Fb are the dominant forces. Since the crustal force is much weaker, 
Tg << Tg, i.e. the evolution of Be is decoupled from that of Be- 

The evolution of core and surface MF and of the spin period for a strong initial MF, 
Beo = lO^"' G, is shown in flgure 5. In this case the strong crustal force prevents the flux 
expulsion from the core much more effectively. The strong and long living MF causes 
an efficient spin-down, up to 30 — 300s, depending on the impurity concentration. The 
residual fleld strength is about 3 ■ 10^ G for Q = 1, and about 3 ■ 10^ G for Q = 0.1. For 
even lower impurity concentrations Tg exceeds the Hubble time. Note that the large initial 
MF enforces Tg = r^. This is a consequence of the strong crustal force, which decelerates 
effectively the expulsion and couples the evolution of the core with that surface MF. 
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Generally, the residual field strength is anti-correlated with the initial field strength and 
positively correlated with Q. Once the balance of crustal and buoyant forces is replaced 
by the balance of vortex acting and buoyant forces on a much lower level, both Be and 
Bf. decay down to a residual value, determined by the final rotational period Py of the 
neutron star. For such neutron stars which do not reach their Pf even after 10^*^ years 
(as for Beo = 10^^ G, Q = 0.01 and Beo = 10^^ G, Q = 1) a residual MF is not attained. 
Among the processes considered here are two dissipative ones: the work done by the drag 
force and the Joule heating produced by field decay in the crust. Assuming that all heat 
produced in that way is irradiated from the surface one can estimate the corresponding 
surface temperature Tg by use of the relation 

Q = A7rR^asBTf , (19) 

where ass is the Stephan-Boltzmann constant, Q = F^-vp + Qjouie, and Qjouie is given by 
equation (12) (see Miralles et al. 1998 for the discussion of the validity of this equation). In 
figure 6 the temporal evolution of Ts is presented for different initial MF strengths and the 
values of the impurity parameter considered above. When the standard cooling scenario 
for a FP-neutron star applies it becomes clear that the contribution of those dissipative 
processes is considerable during the photon cooling era in relatively old (t > 10^ years) 
radio-pulsars. 

4 Discussion 

We considered the effect of the neutron star crust onto the expulsion of a core MF and its 
ohmic decay in the crust. To this aim we solved self-consistently the equation of balance 
of the powers of forces acting on the fluxoids in the core and the rate of change of magnetic 
energy outside the core of the neutron star, assuming a homogeneous MF in the core for 
all the life of the neutron star. 

The evolution of B^, determines the rotational evolution of the neutron star, which in turn 
has considerable effects on the vortex acting force defining the residual field strength. 
It turns out that the characteristic time-scale of the decay of the surface MF r^ increases 
with increasing initial MF strength and decreasing impurity parameter. The amount of 
i?e^decay at t > r^ is correlated with the spin-down rate. A slow spin-down to P/ ~ 1 
s results in a small field decay by less than two orders of magnitude (Pgo = 10^^ G, 
Q = 0.01) while a drastic spin-down to Pf ~ 300 s yields a field decay by about seven 
orders of magnitude (Peo = 10^^ G, Q = 0.1). 

Comparing our results with those obtained for a purely crustal field decay in isolated 
neutron stars (Urpin & Konenkov 1997), we find a considerable deceleration of the decay 
of a field penetrating the entire star. The field evolution for the FP model with standard 
cooling and Q = 0.01 in the case of crustal MF yields an impurity dominated decay with 
Tg ~ 10® years which becomes then nearly exponential. In the present model, only for the 
small initial field of 10^^ G, Tg ~ 10® years, while it is 10^ years for Pgo = 10^^ G and is in 
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the order of the Hubble time scale for B^o = 10^"^ G. Evidently, for a purely crustal MF 
with the inner boundary condition s{Rc, t) = valid for all neutron star's life, no residual 
field can be obtained. The existence of a residual field as well as the strong dependency of 
Tg on BgQ reflects the nonlinear mutual dependency of the field and rotational evolution, 
governed by the balance of forces acting upon the fluxoids. 

Bhattacharya & Datta, 1996, studied the decay of a neutron star MF just expelled from 
the core and deposited in the bottom layers of the crust. They found a rather strong 
decrease of the final MF strength after 10^° years with increasing impurity parameter. 
This result can not be confirmed by our investigation: a larger Q results in a higher 
residual field. This is due to the fact that for a larger Q the crustal force is less strong, 
the time-scale of expulsion is shorter, the surface field decays faster, the spin-down is less 
efficient and the residual field determined by balance of the F„ and Fb is stronger. 
Taking into account that under the assumption of a purely crustal MF a value oi Q = 
0.1. ..0.01 is at least not in disagreement with observations (Urpin & Konenkov 1997), 
these values taken for our model result in a constant B^ for almost all conceivable pulsar 
lifetimes; even for Beo = 10^^ G and Q = 1 a remarkable field decay would start only 
for t > 10^ years. Note that the early (t < 10^ years) cooling-determined crustal field 
evolution is included in our investigations. Only extremely large impurity parameters 
(perhaps Q > 1, which would refiect qualitative deviations from the bcc crystalline struc- 
ture of the crust) would allow for a field decay in t < 10^ years. Thus assuming an initial 
field strength larger than 10^^ G, our study yields a r^ > 10^ years, in agreement with 
the statistical results for isolated radio-pulsars found by Bhattacharya et al. (1992) and 
Hartman et al. (1996). 

It is clearly seen from figure 3, that the main force which is responsible for expulsion of the 
fiux from the core is the buoyancy force. The hypothesis about the so called " spin-down 
induced" expulsion of magnetic fiux (Konar, Bhattacharya 1998 and references therein) 
seems to be adequate only in case of a weak (~ 10^^ G) initial magnetic field. If the 
magnetic field is stronger, say, 10^^ G, vortices will cut through the fiuxoids, whereas the 
latter are anchored in the crust for the time of > 10^ years even in case of high impurity 
concentration [Q = 1). During this time the neutron star spins down by a factor of about 
100, while Be remains almost the same, in contradiction with the core field evolution 
predicted by the spin-down induced mechanism of expulsion. 

In the present paper we study the magneto-rotational evolution of isolated neutron stars 
suffering a spin-down by magneto-dipole radiation. We found that for relatively large 
field strengths, lO^^G < B^q < lO^^G, the residual field -Brcs oc Pj^-, in accordance with 
the result of DCG. We also found that in this range of -Bgo the expulsion times-scale 
is determined mainly by the balance of buoyant and crustal forces. These forces are 
independent of the spin-down rate of the neutron star, i.e. the expulsion time-scale does 
not depend on the specific braking mechanism. Jahan Miri & Bhattacharya, 1994, studied 
the evolution of neutron stars in binaries. They adopted the hypothesis on the spin-down 
induced magnetic fiux expulsion and found B^^s oc P^axi where Pmax is the maximum 
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period reached by the neutron star during the propeller phase. We do not expect that 
our model (if applied to the neutron star evolution in binaries) will confirm this result. 
In the context of the current discussion of magnetars, which are thought to be highly 
magnetized isolated neutron stars spinning down by magneto-dipole braking, the inves- 
tigation of the crustal effect onto the expulsion of initial MFs > lO^^G is an urgent task. 
Our present model is not suitable for that purpose since we consider only ohmic diffusion 
in the crust. In the case of extremely large field strength as expected for magnetars, the 
field evolution in the crust is more complicated, effectively accelerated by a Hall cascade 
(Goldreich & Reisenegger 1992) and/or fracturing of the crust (Thompson & Duncan 
1996). Thus, in order to describe the flux expulsion in magnetars our model has to be 
modifled qualitatively. 

Both the investigation of flux expulsion in magnetars and the effects of accretion on the 
evolution of a MF permeating the whole neutron star, will be considered in forthcoming 
papers. 
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Figure captions 

Fig. 1 Due to the motion of fluxoids from position 1 to position 2 the magnetic field hnes 

are bended in the crust. This causes a change of the magnetic energy and heat release 

within the crust. Thus, a force raises which can counteract the movement of fluxoids. 

Fig. 2 The same as in figure 1 but in terms of the normalized stream function s{r,t). 

The gradient of s near the crust-core interface corresponds, according to equation (9), to 

a generation of a ^-component of the crustal magnetic field. 

Fig. 3 The temporal evolution of magnetic field, forces, velocities of both kinds of tubes 

and spin period of the neutron star for B^o = 10^^ G. Left, middle and right columns 

correspond to the different values of the impurity parameter: Q = 1, 0.1, 0.01, respectively. 

Fig. 4 The same as in figure 3 but for B^q = 10^^ G. 

Fig. 5 The temporal evolution of the magnetic field and spin period for B^o = 10^^ G. 

The numbers at the curves correspond to the different values of the impurity parameter 

Q. 

Fig. 6 The temporal evolution of the surface temperature for different values of B^q and 

Q. Curve 1: no additional heating, curve 2: B^o = 10^^ G, Q = 0.1; curve 3: B^o = 10^^ 

G, Q = 0.01; curve 4: B^o = 10^^ Q, Q = 0.01; curve 5: B^o = 10^^ Q, Q = 0.1; curve 6: 

Beo = 1013 G, g = 1. 
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